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Abstract 

Quality by design (QbD) concept is a paradigm for the improvement of botanical injection quality control. In this work, 
water precipitation process for the manufacturing of Xueshuantong injection, a botanical injection made from Notoginseng 
Radix et Rhizoma, was optimized using a design space approach as a sample. Saponin recovery and total saponin purity 
(TSP) in supernatant were identified as the critical quality attributes (CQAs) of water precipitation using a risk assessment for 
all the processes of Xueshuantong injection. An Ishikawa diagram and experiments of fractional factorial design were 
applied to determine critical process parameters (CPPs). Dry matter content of concentrated extract (DMCC), amount of 
water added (AWA), and stirring speed (SS) were identified as CPPs. Box-Behnken designed experiments were carried out to 
develop models between CPPs and process CQAs. Determination coefficients were higher than 0.86 for all the models. High 
TSP in supernatant can be obtained when DMCC is low and SS is high. Saponin recoveries decreased as DMCC increased. 
Incomplete collection of supernatant was the main reason for the loss of saponins. Design space was calculated using a 
Monte-Carlo simulation method with acceptable probability of 0.90. Recommended normal operation region are located in 
DMCC of 0.38-0.41 g/g, AWA of 3.7-4.9 g/g, and SS of 280-350 rpm, with a probability more than 0.919 to attain CQA 
criteria. Verification experiment results showed that operating DMCC, SS, and AWA within design space can attain CQA 
criteria with high probability. 
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Introduction 

Quality by Design (QbD) is a systematic approach based on 
knowledge management and risk management [1-3]. It has 
become a paradigm for pharmaceutical industry to establish 
analytical methods [4-6], develop new manufacturing processes 
[7-12], and optimize existing manufacturing processes [13]. 
Manufacturing process optimization under the framework of 
QbD will benefit the target patient population, the pharmaceutical 
industry, and the regulation agency. Implementation of QbD 
concept usually contains several steps, such as critical quality 
attribute (CQA) definition, risk assessment, critical process 
parameter (CPP) determination, design space development, 
control strategy design, and continual improvement in drug 
lifecycle [13-15]. 

To keep batch-to-batch consistency, design space development 
is very important in the implementation of QbD concept. 
Parameter variations within design space are not considered to 
affect product quality. Quantitative models between process 
parameters and process CQAs should be established to develop 
design space. CPPs can be determined using risk assessment or 



design of experiments. Plackett-Burman design or fractional 
factorial design are usually applied to select CPPs [16]. In order 
to obtain design space, response surface methodology, such as 
central composite design [17] or Box-Behnken design [18], is 
widely applied. Recently, Rozet et al. suggested that design space 
is "a multivariate domain of input factors ensuring that critically 
chosen responses are included within predefined limits with an 
acceptable level of probability" [16]. The Monte-Carlo method is 
adopted by several researchers to calculate the probability in the 
development of analytical methods [19,20]. 

Botanical drugs have gained increasing popularity in recent 
years. The compositions of most botanical drugs are very 
complicated. There is usually more than one distinct active 
ingredient in a botanical drug. Accordingly, the quality control of 
botanical drugs is difficult. Recently, Zhang et al. applied the QbD 
concept to develop botanical drug manufacturing processes [15]. 

In this work, water precipitation process in the manufacturing of 
Xueshuantong powder was optimized in the framework of QbD 
concept as a sample. Xueshuantong powder was applied to treat 
retinal vein occlusion and cerebrovascular disease in China [21]. It 
is made from Notoginseng Radix el Rhizoma (Sanqi) through a 
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Figure 1. The manufacturing processes of Xueshuantong powder. 
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series of unit operations, including extraction with a mixed 
ethanol-water solution, concentration, water precipitation, decol- 
orization, column chromatography, and freeze drying, as shown in 
Figure 1 [22]. The active pharmaceutical ingredients (APIs) of 
Xueshuantong powder are saponins, including Notoginsenoside 
Ri, ginsenoside Rgi, ginsenoside Rbi, ginsenoside Rd, and 
ginsenoside Re. 

The mixed solvent of ethanol and water is widely applied to 
extract active ingredients from medicinal herbs [23-26]. Water 
precipitation is a precipitation process usually used in the 
treatment of plant extracts extracted using ethanol solution. It 
can be realized by simply adding water into the concentrated 



Table 1. Coded and uncoded values for the factors. 



extracts of plants. The purpose of this process is to remove 
impurities with weak polarity. In the production of Xueshuantong 
injection, the loss of saponins are observed in water precipitation. 
Therefore water precipitation process affects both drug efficacy 
and drug safety. However, there are very few publications on 
water precipitation. 

In this work, water precipitation was optimized using a design 
space approach. Process CQAs of water precipitation were 
obtained through risk assessment. Process parameters are also 
discussed in risk assessment. Fractional factorial designed exper- 
iments were used to determine CPPs of water precipitation. Design 
space was calculated based on the Box-Behnken designed 







Parameters 


Symbols 


Coded values 








-1 


1 


DMCC (g/g) 


A 


0.34 


0.55 


AWA (g/g) 


B 


2.5 


5.5 


Stirring speed (rpm) 


C 


200 


400 


Refrigeration temperature ( C C) 


D 


5.0 


30.0 


Refrigeration time (h) 


E 


12 


48 


Flowrate (ml/min) 


F 


3.0 


7.0 
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Table 4. Conditions and results of verification experiments. 



No. 




VI 


V2 


Concentrated extract amount (g) 




3.0 


3.0 


DMCC (g/g) 




0.52 


0.41 


AWA (g/g) 




2.5 


3.7 


SS (rpm) 




400 


320 


Probability 




0.00 


0.95 


Within design space 




No 


Yes 


Notoginsenoside Rt recovery {%) 


Experimental value 


84.1 ±1.1 


89.2±0.5 




Predicted value 


82.3 


88.6 


Ginsenoside Fig, recovery (%) 


Experimental value 


84.6 ±1.0 


87.1 ±0.2 




Predicted value 


79.3 


87.7 


Ginsenoside Re recovery (%) 


Experimental value 


83.9±1.9 


87.9±0.6 




Predicted value 


79.9 


88.6 


Ginsenoside Rb n recovery (%) 


Experimental value 


89.1 ±1.1 


91.7±0.2 




Predicted value 


83.1 


90.3 


Ginsenoside Rd recovery (%) 


Experimental value 


90.5 ±1.4 


92.3±0.5 




Predicted value 


86.4 


90.0 


Total saponin purity in supernatant (%) 


Experimental value 


40.4±0.2 


40.3±0.5 




Predicted value 


40.3 


41.5 
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experimental results with a Monte-Carlo simulation method. 
Finally, design space was verified. 

Experimental 

Materials and Chemicals 

Sanqi was collected from Wenshan of Yunnan Province (China) 
and purchased from Zhejiang Chinese Medical University 
Medical Pieces Co. Ltd. (Hangzhou, China). No specific 
permissions were required for the described field studies. The 
locations are neither privately owned nor protected by the Chinese 
government. No endangered or protected species were sampled. 

Standard substances of notoginsenoside R^ ginsenoside Rg 1; 
ginsenoside Re, ginsenoside Rbi, and ginsenoside Rd were 
purchased from Shanghai Winherb Pharmaceutical Technology 
Development. Co., Ltd. (Shanghai, China). Acetonitrile (HPLC 
grade) and methanol (HPLC grade) were obtained from Merck 
(Darmstadt, Germany). The formic acid (HPLC grade) was 
purchased from Tedia (Darmstadt, Germany). Deionized water 
was produced using a Milli-Q, academic water purification system 
(Milford, MA, USA). 

Procedures 

Sanqi was crushed and sieved. After that, a mixed solvent of 
60% ethanol and 40% water (v/v) was used to extract Sanqi with 
water bath (W501, Shanghai Shen Sheng Biotechnology Co., Ltd) 
operated at 95°C for three times. The obtained extracts were 
mixed and concentrated. The dry matter content of concentrated 
extract (DMCC) was 509. 1 mg/ g. Concentrated extracts with 
lower DMCC were obtained by dilution with water. 

Water was added to concentrated extract in a conical flask 
under stirring with a roller pump (BT300-2J, Baoding Longer 
Precision Pump Co., Ltd.) at a desired flowrate. After adding 
water, the stirring was kept for 10 min. The flask then was 
refrigerated in a low-temperature thermostat bath (THD-1008W, 



Ningbo Tianheng Instrument Factory) for the desired time at the 
desired temperature. Finally, the supernatant was collected and 
weighed. Saponin contents and dry matter content of the 
supernatant then were determined. 

Experimental design 

Six parameters were investigated, including DMCC, amount of 
water added (AWA), stirring speed (SS), refrigeration temperature, 
refrigeration time, and the flowrate of water addition. The coded 
and uncoded values of parameters are listed in Table 1 . Fractional 
factorial designed experiments were applied to select CPPs, as seen 
in Table 2. Box-Behnken designed experiments were applied to 
obtain the quantitative models between CQAs and CPPs, as seen 
in Table 3. After the development of design space, verification 
experiments with conditions listed in Table 4 were carried out. 
Verification experiments were repeated three times. 

Analytical methods 

The quantitative analysis of saponins was performed using a 
UPLC method. A Waters Acquity UPLC system (Waters, Milford, 
MA) coupled with a Waters Acquity CSH C18 column 
(50 mmx2.1 mm, 1.7 um) was used. The mobile phase consisted 
of 0.01% (v/v) formic acid in water (FA) and 0.01% (v/v) formic 
acid in acetonitrile (FB). The separation was achieved using a 
gradient elution program as follow: 19-20% FB at 0-6 min, 20- 
31% FB at 6-8.5 min, 31-33% FB at 8.5-11 min, 33-90% FB at 
11-17 min, and 90-90% FB at 17-19 min. The solvent flow rate 
was 0.35 ml/min, and the column temperature was kept at 45°C. 
The detection wavelength was 203 nm. Dry matter content was 
determined gravimetrically using a precision electronic balance 
(AB204-N, Metder Toledo Shanghai Co., Ltd). Before weighed, 
samples were dried at 105°C in an oven (DZF-6050, Shanghai 
Jing Hong Laboratory Instrument Co., Ltd.) for 3 h and then kept 
in a desiccator for 0.5 h. 
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Data processing 

The recovery of saponins (SR) are defined as following 
equation. 



SR M SUP -Q >SUP xl0Q% 
Mce'C^ce 

where M and C refer to the mass and concentration, respectively; 
subscript SUP and CE are supernatant and concentrated extract, 
respectively; subscript i (i= 1 to 5) represents notoginsenoside R 1; 
ginsenoside Rgi, ginsenoside Re, ginsenoside Rbi, and ginsenoside 
Rd, respectively. The content of total saponin (Cxs) in supernatant 
or concentrated extract was calculated using Equation 2. 



Cts=J2Q (2) 

/=i 

The purity of total saponin (TSP) was calculated using Equation 3. 



rSP=-^xlOO% (3) 
DM 

where DM is the dry matter content of supernatant or 
concentrated extract. 

Design Expert V8. 0.6.1 (State-Ease Inc., MN) was used to 
analyze the results of fractional factorial design and Box-Behnken 
design. In the selection of CPPs, only the main effects of factors 
were considered, as shown in Equation 4. 



Y = ao + aiA -\- ajB -\- cit,C + a$D + a^E + a&F (4) 

where Y is the response; A, B, C, D, E, and F represent 
parameters; a 0 is a constant; a t to a 6 are regression coefficients. 
The significance level was set to 0.05. Equation 5 was used to 
model the results of Box-Behnken experiments. 

Y=b 0 +biA+b 2 B+b 3 C+b 4 AB+b 5 AC+ 

(5) 

b 6 BC + b 7 A 2 + h B 2 + b 9 C 2 

where b 0 is a constant, and b! to b g are regression coefficients. All 
the parameters were coded before modeling. 

The design space was calculated using a Monte-Carlo method 
with a self-written program of Madab (R20 1 Ob, Version 7.11, 
MathWorks, USA). Considering the sampling error, saponin 
contents and dry matter contents in supernatant were assumed to 
follow a normal distribution. To generate random data of saponin 
contents and dry matter contents, the average values and standard 
deviations were required. The experimental values of saponin 
contents and dry matter contents were used as the average values 
in simulation. The relative standard deviations (RSD) of saponin 
contents and dry matter contents of supernatants were considered 
to be the same with the RSD values of the center point in 
experimental design. Then the standard deviations can be 
calculated according to average values and RSD values. Coded 
values were used in the calculation. Therefore the ranges of 
DMCC, AWA, and SS were within -1 and 1. The calculation 
step sizes for DMCC, AWA, and SS were 0.02, 0.02, and 0.02, 
respectively. Saponin recovery and TSP in supernatant were 
calculated in each simulation. Simulation was repeated 50000 
times to calculate reliable probability values. The acceptable level 
of probability for design space was set as 0.90. 
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Table 6. Criteria of process CQAs. 





CQAs 


Lower control limit (%) 


Upper control limit (%) 


Notoginsenoside recovery 


85 


92 


Ginsenoside Rg 7 recovery 


85 


91 


Ginsenoside Re recovery 


85 


94 


Ginsenoside Rb n recovery 


85 


95 


Ginsenoside Rd recovery 


85 


99 


TSP in supernatant 


40 


43 
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Supernatant remaining in the pores or surface of precipitation 
(SRP) was not collected. The mass of SRP (M SRP ) was calculated 
using Equation 6. 



M S RP = M CE + M Wa ,er - M SU P ~ 

{DMCC x M CE - C dm ,sup x M sup ) 



(6) 



The ratio of saponins in SRP (SRSRP) was calculated using 
Equation 7. 



SM f, = ^ X ^ W xl00% 
M C£ x C, .c£ 



(7) 



Results and Discussion 

Process CQA determination 

There are many CQAs for Xueshuantong powder, including 
color, water content, pH value, saponin purity in dry powder, 
fingerprint similarity, abnormal toxicity, residue on ignition, 




( ^Environment^ ) Q Equipment ^ ^ Refrigeration^ ) 



Figure 2. Ishikawa diagram of water precipitation process. 

doi:1 0.1 371 /journal.pone.01 04493.g002 
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Figure 3. Pareto chart of parameters. 
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(a) Notoginsenoside R] (b) Ginsenoside Rg] 




(c) Ginsenoside Re (d) Ginsenoside Rbi 




(e) Ginsenoside Rd 

Figure 4. Effects of DMCC and AWA on saponin recoveries. 
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Figure 5. Effects of DMCC, and SS on TSP in supernatant. 

doi:10.1371/journal.pone.0104493.g005 



insoluble particles, bacterial content, heavy metals, and harmful 
elements. Heavy metals and harmful elements are generally 
considered to be affected by raw material quality. Bacterial 
content was controlled by producing Xueshuantong powder in 
sterile conditions. Process influences on other drug CQAs are 
identified using risk assessment based on experiences, as seen in 
Table 5. As one of the unit operations in the manufacturing of 
Xueshuantong powder, water precipitation is considered to affect 
saponin purity in dry powder and fingerprint similarity. Higher 
active ingredient recovery and higher TSP in supernatant both 
help to realize higher saponin purity in dry powder. Higher TSP in 
supernatant also means higher removal of impurities, which affects 
fingerprint similarity. Therefore saponin recovery and TSP in 
supernatant were identified as the process CQAs in this work. 
Their criteria are listed in Table 6. 

Selection of CPPs 

An Ishikawa diagram analysis was performed to find out process 
parameters that affect the CQAs of water precipitation, as seen in 
Figure 2. The risk severity of four parameters of concentrated 
extract temperature, concentrated extract pH value, environmen- 
tal temperature, and water addition style is low. Parameter fail 
probability of water addition style, and style of mixing is very low. 
The detectability for concentrated extract temperature, and 
environmental temperature is high. Therefore, the risk caused 
by these five parameters are low. 

The other six parameters of AWA, SS, DMCC, refrigeration 
temperature, refrigeration time, and water addition flow rate were 



investigated using fractional factorial design. The results of 
fractional factorial design experiments are listed in Table 2. The 
Pareto chart for saponin recovery and TSP in supernatant are 
shown in Figure 3. In Figure 3(a) and 3(b), DMCC and AWA 
remarkably affected the recoveries of notoginsenoside Rj and 
ginsenoside Rgi. Ginsenoside Re recovery was mainly affected by 
AWA, as seen in Figure 3(c). In Figure 3(f), SS significandy 
affected TSP in supernatant. The other three factors, including 
refrigeration temperature, refrigeration time, and the flow rate of 
water were insignificant on process CQAs. Therefore DMCC, 
AWA, and SS were selected as CPPs. 

The effects of CPPs on process CQAs 

The experimental results of Box-Behnken design are displayed 
in Table 3. The recoveries of all the five saponins are higher than 
78%, which means most active ingredients can be recovered in 
water precipitation. TSP in supernatant varied from 38.8% to 
43.2%. Because TSP in concentrated extract was 36.7%, it can be 
concluded that TSP increases after water precipitation. Response 
surface models were developed to obtain quantitative relationships 
between CQAs and CPPs. The estimated regression coefficients 
are listed in Table 7. Analysis of variance (ANOVA) was 
performed, and p values of parameters are also listed in Table 7. 
For all the models, determination coefficients (R 2 ) are larger than 
0.86, which means most variations of process CQAs can be 
explained by DMCC, AWA, and SS. Models are significant with 
p-values less than 0.05. In Table 7, DMCC are significant for all 
the CQAs. AWA are significant for all the recoveries of saponins. 
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Figure 6. Sum of SR and SRSRP for different saponins. 
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Run 



The quadratic term of DMCC shows significant influences on the 
recovery of notoginsenoside R[ and ginsenoside Rgi- The purity of 
total saponin in supernatant is mainly affected by DMCC and the 
interaction term of DMCC and SS. 

Figure 4 shows the effects of DMCC and AWA on the 
recoveries of saponins. The increase of AWA can reduce saponin 
loss. However, saponin recoveries decrease as DMCC increases. 
The effects of CPPs on TSP in supernatant are shown in Figure 5. 
High TSP in supernatant can be obtained when DMCC is low 
and SS is high. 

Loss mechanism of saponins 

Saponins will degrade under acidic, basic, or thermal conditions 
[27-29]. The solubilities of some saponins were also reported 
[30,31]. It can be concluded that saponin solubility in water is not 
large. Theoretically, the loss of saponins in water precipitation 
may be caused by chemical reactions, precipitation, or the 
incomplete collection of supernatant. The sums of SR and SRSRP 
for each saponin are plotted in Figure 6. Most values of the sums 
of SR and SRSRP are within 100% ±5%. It indicates that the 
incomplete collection of supernatant is the main reason for the loss 
of saponins. The uncollected supernatant mainly remained in the 
pores and surface of precipitation. 



Design space development and verification 

The calculated design space are shown in Figure 7. It can be 
concluded that the design space is an irregular polygon. 
Recommended normal operation region to attain CQA criteria 
are located in DMCC of 0.38-0.41 g/g, AWA of 3.7^4.9 g/g, and 
SS of 280-350 rpm with a probability more than 0.919. 
Verification experiments were carried out. The prediction results 
and experimental results of process CQAs are listed in Table 4. 
The recoveries of notoginsenoside R 1; ginsenoside Rg 1; and 
ginsenoside Re of Experiment VI are lower than their criteria. 
While all the experimental values of Experiment V2 were in the 
ranges listed in Table 6. These results indicate that criteria of 
CQAs can be attained by operating CPPs within design space. 

Conclusions 

In this work, water precipitation process for the manufacturing 
of Xueshuantong injection was optimized using a design space 
approach as a sample. Saponin recovery and the total saponin 
purity in supernatant were identified as the process CQAs of water 
precipitation using risk assessment. An Ishikawa diagram was 
applied to find out potential CPPs. DMCC, SS, and AWA were 
identified as CPPs using a fractional factorial design. Models 
between CPPs and process CQAs were developed with experi- 
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Figure 7. Design space obtained using Monte-Carlo simulations. 
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mental results of Box-Behnken design. Determination coefficients 
were higher than 0.86 for all the models. High TSP in supernatant 
can be obtained when DMCC is low and SS is high. Saponin 
recoveries decrease as DMCC increases. The loss of saponins was 
mainly caused by the incomplete collection of supernatant. Design 
space was calculated with a Monte-Carlo simulation method. In 
the simulation, 0.90 was adopted as the acceptable probability for 
the attainment of CQA criteria. Recommended normal operation 
region are located in DMCC of 0.38-0.41 g/g, AWA of 3.7- 
4.9 g/ g, and SS of 280-350 rpm, with a probability more than 
0.919 to attain CQA criteria. Design space was verified and 
verification experimental results showed that operating DMCC, 



SS, and AWA within design space helps to attain process CQA 
criteria with high probability. The design space applied in this 
work can also be used to optimize other unit operations in 
pharmaceutical engineering. 
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